Cluster interactions for fcc-based structures in the Al-Mg-Si system 
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A class of proposed coherent precipitate structures (Guinier-Preston zones) in the Al-Mg-Si alloy 
are investigated using first-principles density functional theory methods. The cluster expansion 
method is used to extract effective interaction parameters, providing the means for large scale 
energy calculations of alloy structures. The MgiSii Llo structure and structures related to the 
MgsSis P" phase are studied in more detail, and e.g., precipitate/matrix interface energies are 
presented. Using direct first-principles calculations we show that the former phase is dynamically 
unstable and thus must be stabilized by the surrounding Al matrix. Monte Carlo simulations and 
free-energy techniques are used to study the Al rich side of the phase diagram with the current CE 
parameters, and kinetic Monte Carlo simulations are used to study clustering in the disordered phase. 
The implications of our findings are discussed in the framework of classical nucleation theories, and 
we outline possible nucleation mechanisms. 
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I. INTRODUCTION 

The formation of solute-rich precipitates in a super- 
saturated solid solution is a process of great technological 
importance because it gives raise to precipitation hard- 
ening in a wide class of alloys. Among the precipita- 
tion hardened Al alloys, the Al-Mg-Si (6xxx) alloys be- 
long to the most studied with regards to their precipita- 
tion sequences, see Ref. [l|, 0, @ and references therein. 
They often display a complex sequence of metastable 
phases forming under different conditions (time, tem- 
perature and composition). The structures of almost 
all of these phases are known in detail from transmis- 
sion electron microscopy (TEM) studies. However, such 
techniques give no direct information about the ener- 
getics of the phases of interest. First-principles calcu- 
lations can provide vital new information in this respect, 
as has recently been demonstrated in the case of the Al- 
Mg-Si system, see work by Derlet et aZQ and Ravi and 
Wolverton.[3 Thus, by means of solving the Schrodinger 
equation within the density functional theory (DFT), one 
can obtain enthalpies, and study relaxation and binding, 
in phases which are meta-stable and only exist as nano- 
sized particles within the host matrix. A limitation of 
direct first-principles calculations is that, in practice, the 
structure has to be specified in advance. A central and 
unresolved question in the case of the Al-Mg-Si system 
is the structure of the first precipitates that form from 
solid solution, i.e., the Guinier-Preston (GP) zones. 3] 

In order to be able to study coherent precipitate phases 
within the Al matrix, and to model the disordered phase 
and the nucleation process, we have utilized a ternary 
cluster expansion (CE) description of the Al-Mg-Si sys- 
tem. In this model, a set of atoms within a lattice struc- 
ture is denoted by <t, and the total energy of such a con- 
figuration can be written as a sum of "cluster functions" , 
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We derive the cluster interaction parameters Va^s from 
a set of first-principles calculations using an approach 
based on the structure inversion method (SIM).^ Once 
these parameters have been obtained, it is a relatively 
inexpensive operation to sum up the total energy. A 
complication that arises is that the effective cluster in- 
teractions are in principle volume dependent, |7| but since 
the concentration of Mg and Si in typical Al-Mg-Si alloys 
is 1 — 2%, wc have done all calculations at volumes cor- 
responding to the Al equilibrium lattice constant. 

We have focused on the two main candidates for early 
precipitate structures, or GP-zones, in the Al-Mg-Si sys- 
tem. These are, the MgiSii LIq structure proposed by 
Matsuda et al,^ here referred to as the Matsuda struc- 
ture, and the so called pre-/3"-phase which was proposed 
by Marioara et al. @ By testing a large number of struc- 
tures related to the pre-/?" structure, we conclude that 
the Matsuda-phase is the energetically most stable pre- 
cipitate phase. We then proceed with studying the sta- 
bility of the Matsuda and pre-/?" structures and find that 
the former structure is dynamically unstable, so that the 
vibrational free energy is undefined. This implies that 
it can only exist within a surrounding stabilizing lat- 
tice, the Al lattice in this case. That this system ex- 
hibits instabilities is not surprising since fee Si is dynam- 
ically unstable. [Toll For certain local concentrations of Si 
we therefore expect the lattice to be significantly soft- 
ened, which lowers the Gibbs free energy, and may lead 
to structural transformations. 

The CE method also allows for large-scale Monte Carlo 
(MC) simulations to be performed, either for finding 
ground-state structures (using, e.g., simulated annealing) 



or to study finite-temperature thermodynamic proper- 
ties of an alloy. Our MC calculations using the present 
CE parameters show that decomposition typically occurs 
in the form of spherical precipitates with the Matsuda 
structure. MC simulations in combination with ther- 
modynamic integration are used to calculate the solvus 
phase boundary in the Al rich end of the phase dia- 
gram. We also use kinetic Monte Carlo (kMC) simu- 
lations to study clustering in the disordered phase. Such 
studies are of current interest because modern experi- 
mental techniques, i.e., high- resolution transmision elec- 
tron microscopy and 3D atom probe tomography, are ca- 
pable of studying small, coherent clusters at an atomic 
scale. l^jlUi Therefore, increased understanding of atomic 
level processes that preceed nucleation can be expected 
and a much more detailed comparison with modelling 
and simulation is possible., llj 

The paper is disposed as follows. In Sec. II we describe 
the two GP-zone models considered in this work. In Sec. 
Ill we review the theory of ternary CE, and we also detail 
the first-principles calculations that form the basis for the 
present work. The results of the fitting procedure, and 
of the CE calculations for different GP-zone structures 
are presented in section IV. The calculation of the solvus 
boundary for MgiSii and simulations of the early stage 
of clustering are presented in the same section. In Sec. 
V we discuss the present results in the light of what is 
experimentally known about GP-zone formation in Al- 
Mg-Si, and conclusions follow in Sec. VI. 



II. GP-ZONE MODELS IN THE AL-MG-SI 
ALLOYS 



Currently there are two generally accepted models of 
the early precipitate phases, or GP-zones, in the Al-Mg- 
Si alloys. The first structure, an MgiSii Llo phase (see 
FigD^) was observed by Matsuda and co-workers p| in 
samples with a solute concentration of 1.55 at.% and with 
a Mg:Si ratio of 2:1. The samples had been aged and 
then heat-treated at 343 K. This phase will be referred 
to as the Matsuda phase. The other structure, the pre-/3" 
in Fig. ^ was proposed by Marioara et al. [3 based on 
TEM observations of samples that had been homogenized 
at 843 K and heat treated at 450 K. The Mg:Si ratio in 
these samples was 5:6. 

The /?" structure is the most efficient hardening phase 
in the 6xxx alloys. The fact that GP-zones have been 
observed to have the same symmetry as the /3" phase, 
motivated the structural model of GP-zones shown in 
Fig. ^. Q Its structural relationship with the j3" phase 
can be understood in the following way. If the (0, 0, 0) 
Mg atom in Fig. QJi is displaced l/2ao in the normal 
(z) direction, one obtains a structural model of the (3" 
phase. The latter phase is in turn needle-shaped, and 
often observed to be coherent in the (a; y) plane and 
semi-coherent in the z (needle) direction. 
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FIG. 1: Previously proposed GP-zone structures, where large 
(small) circles indicate Mg (Si) atoms. Dark atoms are dis- 
placed by 1/2 ao in the z-direction, with respect to the white 
atoms, a) the pre-/3" structure with composition MgsSis, the 
11 atom primitive unit cell is indicated, b) the Matsuda Llo 
structure with composition MgiSii. 



III. METHODS 



Cluster expansion and the structure inversion 
method 



The general theory for multicomponent CE was de- 
veloped by Sanches, Ducastelle and Gratias.flJl In the 
present work we use the following notations. A lattice 
site i can be occupied by an atom of type A^ B, or C, 
indicated by an occupation variable, ai, taking the value 
— 1, 0, or -f-l, respectively. Hence, given an underlying 
lattice, the alloy structure can be defined by a set of oc- 
cupation variables, or = {(Ti}. Linearly independent basis 
vectors describing the occupation are then 
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A cluster a is defined as a set of lattice sites 
(pi,P2, ••■,Pn)- Given a cluster a, a cluster function is 
defined as the the product of the basis vectors. 
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where the sequence s = (si, S2, ..., s„) is called "decora- 
tion", which in the ternary case will consist of l:s and 
2:s. The total energy is given by (Q, and by considering 
the average of cluster functions for symmetry-equivalent 
cluster-decoration pairs, ^a,s{o'), the total energy can be 
re- written as 



E{<t) = Eo ■ 
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where Eq is a reference energy and $q^s(c) is an average, 
defined as 



^a,.(o-) = — ^$a',,'(cr) 



(5) 
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The sum goes over all {a',s') pairs related to (a, s) 
through lattice symmetry operations, rua^s is the num- 
ber of such cluster functions per lattice site, and N is the 
total number of lattice sites. 

In the structure inversion method one calculates the 
VqIS by minimizing 
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where {cr} denotes a list of n input structures and 
_Edft(c) their corresponding total energies. The choice 
of input structures is arbitrary as long as they contain in- 
formation about all cluster functions, and an advantage 
with the SIM is that the cr:s can consist of small unit 
cells in a super-cell approach, for which the total energy, 
may be efficiently calculated using DFT methods. 

In order to select cluster functions to be used in the 
fitting, we used cross-validation, 13j with the intention 
to obtain an optimal fit, while at the same time avoid 
over-fitting. It is applied to CE by defining a hierarchy 
of clusters, were a cluster is included only if all its sub- 
clusters are included, and short range clusters are priori- 
tized over long-range clusters. Then, by choosing clusters 
in this way, the aim is to minimize the cross-validation 
score 



(cv)' 



^E(^'-%))^ 



(7) 



Here, _B(i) denotes the predicted energy of structure i, us- 
ing parameters V^.s obtained without structure i included 
in the fit. Thus, the cv-score measures the ability of the 
CE, with a given set of cluster functions, to predict the 
energy of new structures. We note that all n structures 
are included in the fit in Eq. @. 



B. Electronic structure calculations 

We have used the Vienna ab-initio simulation pack- 
age (VASP),|lJ 3. plane- wave pseudo potential imple- 
mentation of the DFT method, in the present electronic 
structure calculations. For the exchange-correlation ap- 
proximation we used the Perdew-Wang implementation 
of the generalized gradient approximation (GGA).|il5ij 
The built-in ultra-soft pseudo potentials by Vanderbilt 
et al. were used. 16] /c-points were chosen according to 
the Monkhorst-Pack scheme in meshes corresponding to 
12 X 12 X 12 points for the conventional fee unit cell, or 
denser. Tests showed that using a mesh corresponding 
to 16 X 16 X 16 fc-points changed the energy with typi- 
cally 1-2 mcV/atom. The atomic structure was relaxed at 



constant volume using a conjugate gradient technique, or 
in some cases a damped molecular dynamics algorithm. 
We found that both methods led to the same relaxed 
state. We note that during the relaxation, the symme- 
try is preserved. This means that if the relaxation starts 
at a saddle point in the potential energy landscape, the 
system may stay at this saddle point due to symmetry 
restrictions. In particular, pure Si in the fee structure 
is unstable with respect to tetragonal shear, lOj but this 
relaxation mode is not taken into account in our calcula- 
tions. 



C. Monte Carlo and kinetic Monte Carlo 
simulations 

The computational study of the thermodynamic prop- 
erties of the Al-Mg-Si system requires to take configu- 
rational space properly into account. To this end we 
use the standard Metropolis 17] Monte Carlo sampling 
technique, generating a sequence of configurations being 
Boltzmann-distributed with respect to the system tem- 
perature and configurational energy. We have held the al- 
loy constituent concentrations unchanged during the sim- 
ulations. Hence, as the trial-changes we simply choose a 
random pair of atoms to be exchanged with each-other. 
One MC step is defined as one such trial-echange per 
atom. A technical aspect of our implementation is that 
the cluster functions are defined in the so called canonical 
basis rather than the orthogonal basis (Eq:s Q). That 
means that the cluster functions are products of integers 
rather than real numbers, and their evaluation becomes 
more efficient. 

In order to simulate kinetic processes (e.g., nucleation) 
we have utilized a kinetic Monte Carlo (kMC) algorithm. 
The important additional constraint added to a kMC al- 
gorithm, as compared to the ordinary MC dito, is that 
the (trial-) changes made to the system are selected in 
a way correspondig to the effect of some physical pro- 
cess occurring. The physical process being most impor- 
tant during aging or phase transformations in our Al- 
Mg-Si alloy system, is the mass transport via atomic dif- 
fusion. The atomic diffusion process was modelled by 
attempted interchange of two randomly choosen nearest 
neighbor atoms, and accepting the interchange according 
to the Metropolis algorithm. This model is also known 
as Kawasaki-dvnamics. [iS) 

To provide an estimate of the time-scale of the sim- 
ulated diffusion process, we use an experimental value 
of the aluminum self diffusion coefficient 19] as input. 
The experimental input which implicitly contain infor- 
mation about diffusion hop-rates and vacancy concentra- 
tions, will relate the Monte Carlo steps to physical time 
as 
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where r is the physical time scale of giving all the atoms 



in the system a chance to diffuse, oq being the Al lattice 
contant, and -Dcxp(T') is the experimentaUy determined 
diffusion coefficient. The factor Sv accounts for the su- 
persaturation of vacancies that is present, e.g., after a 
quench. 

We used the kMC method to study the early stage of 
clustering in Al-Mg-Si. Our definition of a cluster is then 
that each solute atom (Mg or Si) belonging to a given 
cluster is nearest neighbour to at least one more solute 
atom in the same cluster. 



D. Free energy integration 

The phase diagram, and the underlying free-energy dif- 
ferences between various phases determine the driving 
force, e.g., for precipitate formation. In order to calcu- 
late the phase boundary between disoerdered and a two- 
phase region, one must know the free energy of the dis- 
ordered phase. In the present work, it was calculated in 
the following way. The energy as a function of tempera- 
ture was determined for a set of compositions in standard 
MC simulations. This was done by slowly decreasing the 
temperature and the system was then kept in the disor- 
derd state. A corresponding simulation where the tem- 
perature is increased from a low temperature shows a 
clear hysteresis and the point where the two curves meet 
can be taken as a rough estimate of the phase boundary. 
The entropy in the disordered state was determined using 
thermodynamic integration of the relation 



dS_ 



1 dH 
T'dT' 



(9) 



The remaining constant of intgration was taken from the 
expression for the entropy of a ternary system in the ideal 
solution model 



Si 



ideal 



-fcB(ci In(ci) + C2 ln(c2) + cg ln(c3)) (10) 

where ci, C2 and C3 are the respective atomic concen- 
trations. This expression is which is valid at sufRciently 
high temperature. In the present work, we integrated 
from T = 3000 K. 



IV. RESULTS 

A. Fitting of the CE parameters 

To choose input structures for the SIM, we started with 
8 binary structures and 34 ternary structures in Refs. 1201 
andl2ll respectively. Those structures have been shown 
to be the only distinct structures that can be obtained if 
one assumes only up to next-nearest neighbor pair inter- 
actions. (In the ternary case, this has not been strictly 
proven, but it seems plausible.) Taking into account the 
structures obtained when the occupation of A, B or C 
atoms in the 'generic' structures are permuted, one ends 
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FIG. 2: Convergence of the cv-score (circles) and the mean- 
square deviation between predicted and calculated energies 
of the input structures (squares), as a function of included 
clusters; pairs (P), three-point clusters (T), and four-point 
clusters (Q). It is understood that, e.g., P4 means Pl-4. See 
further Table 1. 



up with 148 symmetrically distinct structures. In addi- 
tion, we used some low-concentration structures, e.g., an 
isolated impurity in a 32 LP cell, see Appendix A. Fur- 
ther, we excluded some structures, e.g., pure Mg, that 
had a large static (positive or negative) internal pressure. 
In total, we ended up with 154 structures. 

These structures were used as input to the fitting of the 
cluster interaction parameters in Eq. @ . When selecting 
clusters, we minimized the cross-validation score given by 
Eq. Q. In the procedure we included pairs, triplets and 
quadruplets in a hierarchical order, see Fig. [3 and Table 
U for an explanation of the clusters used. We found that 
a good fit was obtained if we included pairs up to fourth 
nearest neighbor, and the triplets Tl, T2 and T4 (see Ta- 
bleOJ. This gives a total of 27 parameters Va^s that were 
minimized. The resulting cv-score is cv = 1 A meV/atom. 



B. Static energies of suggested GP-zone structures 

We have investigated the energetics of two models of 
the GP-zones in Al-Mg-Si; the Matsuda phase and the 
pre-/3" phase (see Sec. II and Fig.^l. 

Because the latter structure has not actually been 
solved, we proceed along the lines in Ref. 5; starting 
with the 11 atom unit cell in Fig. ^, we exchanged 
Mg/Si atoms with Al atoms in all 2^^ = 2048 possible 
ways. For each structure (of which some will be identi- 
cal by symmetry) we calculated the total energy using 
Eq. 1^. In Fig. U] we have plotted, for each composi- 
tion CMg/(cMg + csi), the energy of the most stable struc- 
ture. An even less restrictive procedure is to generate all 
3^^ ~ YIIX'^ structures that arise if each lattice point 
is occupied by a Al, Mg or Si atom. The most stable 
structures in this case are essentially the same as those 
in Fig. 

Ravi and Wolvertonj^ did a set of calculations for pre- 
/3" structures, using first-principles methods. We have 
used the same convention as they used in plotting the 
energies, i.e., we plot the formation energies per solute 
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FIG. 3: Total energy per solute atom of different suggested 
GP-zone structures, along with some of the structures used in 
the fitting procedure. The results are given with reference to 
the energies of the pure elements, in the fee structure, at zero 
pressure. The structures marked with diamonds are denoted 
in accordance with Ref. Hu 



atom with respect to the pure constituents in the fee 
structure. Therefore, the the two sets of energies can be 
directly compared. 

Fig. 121 shows that structures with an Mg:Si composi- 
tion close to 1:1 are energetically favoured. This is in line 
with most experimental data on clusters and early pre- 
cipitates in Al-Mg-Si alloys. 8j For this composition, the 
Matsuda phase is lowest in energy, as found also in Ref.|5i. 
In order to picture how "pre-/3"-like" structures may be 
build up, we also extracted the energetically most stable 
structures contained in the pre-/3" cell as a function of 
solute concentration Csoi = cug + csi — 1/11, ...,11/11. 
Selecting the structures from the 2^^ or the 3"'^^ possible 
structures resulted in the same ground state structures. 
The energies are plotted in Fig. 21 together with directly 
calculated energies for the same structures. Thus, by 
comparing the two sets of energies, one gets an estimate 
of the predictive power of the CE method used. We note 
the the relaxation of the structure at Csoi = 9/11 did not 
finish during the simulation. This could indicate that the 
structure is unstable. Some representative lowest-energy 
structures in Fig. 0] are shown in Appendix B. 

The interface energies between secondary phase par- 
ticles and a parent phase are important parameters in 
various theories of nucleation and phase stability that 
are usually hard to estimate. We calculated the pre- 
f3" /A\ and Matsuda/ Al interface energies (at T = 0) us- 
ing Eq. ((U, see TableHIl The finding that the pre-/3"/Al 
(3 1 0) interface is negative was confirmed in a direct 
DFT calculation in a 44 atom super-cell. The result in- 
dicates that this phase may contain Al. 

The energy E{(t) of Eq. Q is the static energy at con- 
stant volume. To this comes a lowering of elastic energy, 
^^ciast ^ j^g |.Q volume and cell-shape relaxation of these 
phases. For the structures that were used in the fitting 
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FIG. 4: Predicted ground state structures in the 11 atom unit 
cell of the pre-/3" structure (circles) as a function of solute 
concentration. Diamonds indicate the total energies obtained 
in direct DFT calculations for the same structures. 



of the CE parameters, this term can be estimated in the 
following way. Assuming that the phase B relaxes centro- 
symmetrically within an A matrix, which is assumed to 
be elastically isotropic, ^E'^^^^'^ is given bv|2^ 



^^olast ^ l^^j^KsliiKB + ^[iA 



)f^VB 



(11) 



Here, [ia is the shear modulus of the A phase and Kb is 
the bulk modulus of the B phase, where in this estimate 
we used iiTsi^fcc = ^ai- vb is the volume of B phase and e 
is the initial linear misfit of the B phase relative to the A 
phase, e is in turn calculated from the internal pressure 
■p of the different structures using Vb = V^/(l — p/K). 
The results for some structures in Fig. ^ are given in 
Table lini 



C. Dynamical stability of the Matsuda and pre-/?" 
phases 

The dynamic stability of the Matsuda phase and the 
pre-/3" phase was tested by displacing the atoms from 
their relaxed positions in the internally relaxed struc- 
tures. For the Matsuda phase, starting with a 32 atom 
cell, we displaced the atoms in random directions 1/50 
of a nearest-neighbor distance. A subsequent relaxation 
showed that the forces were non-conservative (i.e., di- 
rected away from the initial positions) for the Si atoms 
within the Si plane. 

A similar analysis of the pre-/?" structure shows that 
it is most likely metastable. In a separate calculation, 
we displaced the (0, 0, 0) Mg atom in Fig. QJi by a small 
distance in the z-direction. Also in this case, the force on 
the Mg atom was conservative. 

In the case if the Matsuda phase, we note that, accord- 
ing to the findings in Ref. |^ the structure exists in thin 
(one or a few atomic) (0 1 0) slabs in the Al lattice (see 
Fig-Ct)- For such morphologies, the Matsuda phase will 
be stabilized by the surrounding Al lattice. 
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FIG. 5: Free energies for the disordered phase, calculated 
by thermodynamic integration from a reference state at 3000 
K. The lines correspond to Csoi ~ 1,2, ...,6 at.% solutes, in 
order of increasing slope. The zero energy correspond to the 
reference state of completely separated Llo and Al phases, 
E = CsoI-Eai + (1 — Csoi)£^Lio. The inset shows one of the 
switching runs, at c^oi ~ 2%. The polynomial fit is included, 
and overlaps with the energy points. 



D. Free energy of precipitation and kinetics of 
clustering 

In order to study precipitation and clustering it is nec- 
essary to know the solubility limit of the alloying ele- 
ments. In this sub-section we first describe the calcu- 
lation of the solvus boundary in the Al rich end of the 
phase diagram, and then we present kMC simulations of 
clustering in the disordered phase. 

We have focused on the subset of the phase diagram 
where the Mg:Si ratio is 1:1 and Csoi = cug + csi varies 
between and 5 %. The energy per atom as a func- 
tion of T was calculated for a range of concentrations, 
while staying in the disorderes phase. In these simu- 
lations, an 80 X 80 X 80 atom system was cooled at a 
rate of 0.02 K/kMC step. Based on a polynomial fit to 
the energy, we then used thermodynamic integration to 
calculate the entropy by use of Eqs. Q and H10() . The 
corresponding free energies are plotted in Fig. |S1 The 
free energies are given with respect to the normal state, 
which corresponds to separated Al and LIq phases at K. 
The solubility limit at each temperature corresponds to 
the lowest lying free-energy curve. We have implicitly 
assumed that the ordered phase, the LIq, is completely 
ordered, i.e., does not contain anti-defects or Al atoms. 
We checked the validity of this assumption, again us- 
ing the regular solution model. It turns out that defects 
(typically 1-2% below 500 K), can be neglected in the 
calculation of the solubility limit at the Al rich end of 
the phase diagram. The resulting solvus-line is plotted 
in Fig. El 

Clustering in Al-Mg-Si was studied by kMC. The 
atoms were initially randomly distributed, correspond- 




FIG. 6: Calculated solvus line dividing the disordered and the 
two-phase regions of the phase diagram for which CMg ~ csi = 
Csoi/2 (squares). The dashed line corresponds to an estimate 
based on the regular solution model. 



ing to a previous homogenisation at a high temperature. 
In Fig. we have displayed the cluster size distributions 
as a function of annealing time for a 2% solution at 400 
K. The time-scale was calculated by assuming the Al self- 
diffusion rate Dai = 1-5 x 10-5exp(-1.33/fcB/400) and 
a 10000-fold increase in the number of vacancies, corre- 
sponding to a homogenisation at 750 K. We believe that 
the latter factor is the dominant source of error in esti- 
mating the time-scale in this type of simulations. It is 
seen that the initial stage of clustering is fast (~ 1 min.). 
In practice, since the kinetics is even faster at higher T, 
part of the initial clustering is expected to take place 
already during the quench. 

Since the system is not within the two-phase region, 
there are no nucleation and growth stages. After the 
initial clustering, the cluster sizes fluctuate in accordance 
within 



Ar<=q(i)/iV = exp(-Ai^(i)/fcBT) 



(12) 



where AF{i) is the excess free energy required to form a 
cluster of size i, and N'^'^ is the number of such clusters 
in a system containing N atoms. AF{i) can be written 
as 



AF{i) = iA/™i + Bi^/^a 



(13) 



where A/voi is the volume free energy per atom of a clus- 
ter, a is the cluster/matrix interface free energy and B 
is a geometrical factor. 

A longer kMC simulation (6 x 10^ kMC steps for 
80 X 80 X 80 atoms) showed that the cluster sizes were 
distributed in accordance with Eq. (|12l) for clusters from 
« = 8 up to i = 16, with an effective volume energy 
close to and an interface energy cr = 3 meV/A , i.e., 
considerably lower than the mean interface energy of the 
Matsuda phase at K (Table |Il|). 




10 100 1000 10000 

kMC steps 

FIG. 7: Log-log plot of the number of clusters of a given 
size i, calculated in a kinetic Monte Carlo simulation of an 
80 X 80 X 80 atom system. The horizontal lines are the mean 
number of clusters over a total of 6 x 10^ kMC steps. kMC 
steps have been related to physical time as described in the 
text. 



V. DISCUSSION 

The main purpose of this paper was to apply the CE- 
SIM method to model clusters and GP-zones in the Al- 
Mg-Si system. The interaction parameters Va.s were fit- 
ted to reproduce the static energy E of fee based struc- 
tures in the low solute-concentration regime. By use of 
cross-validation and by selecting clusters of increasing 
range, we ensure that an 'optimal' set of clusters is se- 
lected in the fitting procedure. Explicit calculations on 10 
predicted ground-state structures derived from the prc- 
/3" structure reveal that the method holds its promise; 
the mean square deviation between predicted and di- 
rectly calculated energies is 7.25 meV/atom, to be com- 
pared with the cv-score cv = 7.1 meV/atom in the fitting 
procedure. 

However, the energy of interest for a given cr is the 
Gibbs free energy G = E — TS + pV . The pV term can 
be handled within the CE formalism, 7] but it is much 
more difficult to treat the TS term in a completely gen- 
eral manner. Our explicit calculations for the Matsuda 
phase show that it is dynamically unstable, i.e., the TS 
term is undefined. In general, one expects that when 
the local configuration approaches an unstable structure, 
the vibrational frequencies are lowered (phonon soften- 
ing), so that the structure is stabilized. Thus, this is 
an effect that we expect to be important in the early 
stages of precipitation in Al-Mg-Si. We note that in Al- 
Sc, the inclusion of the TS term changes the phase dia- 
gram significantly. [2 4] This is related to phonon softening 
but not to instabilities. Instabilities in substitutional al- 
loys have been studied previously, e.g. in W-Re.!2a| 

We now discuss the implications of our findings in re- 
lation to what is known from experiments on the nucle- 
ation of GP-zones and the /9"-phase. In our simulations, 
we have only found one stable precipitate structure; the 
MgiSii LIq structure. It forms in approximately spher- 



ical precipitates below T k, 400 K for the composition 
CMg + csi = 2 %. Matsuda et al. report the formation 
of MgiSii Llo precipitates at T = 323 K for an alloy 
containing 1.6 mass % Mg2Si. These precipitates are 
in the form of (0 1 0) planes. The morphology can be 
related to the internal stress of the Matsuda structure, 
which favours platelets in the (0 1 0) plane. Because the 
platelets extend only one or a few planes in the normal 
direction, the Matsuda phase will be stabilized by the 
surrounding Al matrix. 

The pre-/3" and (3" phases are found experimentally 
to form at r = 423... 523 K in alloys containing about 
1.25 at.% Al and Mg. 3] Fig. shows the calculated 
phase boundary for forming Matsuda particles. One finds 
that neither the Matsuda phase nor the pre-/?" phase 
would form with the current atomic interaction parame- 
ters, at the temperatures and solute concentration where 
/3" forms. Elastic relaxations would lower the energies of 
both the impurities and the ordered phases, and are not 
expected to change this conclusion (see further Ref.Q)- It 
is seen in Fig. that the regular solution model describes 
the solvus phase boundary well up to Csoi ~ 4%. For a 
Mg:Si composition of 1:1, the expression for the solvus 
line according to this model is c w 2exp(— Aii^/fce^^), in 
the low concentration limit. AS is the ordering energy 
which can be found in Fig. |21 by comparing the energy 
of, e.g., the Matsuda phase with a line connecting the 
impurity energies of Mg and Si. Thus, one finds that 
the ordering energy required to stabilize a structure at 
the temperature and concentration where GP-zones form 
(w 450 K and 1.5% repectively), is at least 0.2 eV/atom. 
This is larger than any of the ordering energies in the 
current study or in Ref. ^5*. 

Therefore, in order to explain the nucleation mecha- 
nism of GP-zones and the j3" phase, we suggest the fol- 
lowing possibilities. Either, vibrational entropy effects 
or vacancies stabilize an ordered pre-/?" based structure, 
which then transforms into the [3" structure, or the struc- 
tural transition takes place in the disordered phase, in 
a Mg-Si rich cluster with the right local configuration. 
Recent TEM observations clearly indicate that the GP- 
zones are ordered, fcc-based, structures.^ Thus, the for- 
mer explanation seems more likely. In either case, there 
is a substantial gain in energy connected to the transition 
to the /?" structure. 

Finally, we remark that it is likely that quenched-in 
vacancies play a direct role in the nucleation of the j3" 
phase. As pointed out in Ref. ^, vacancies may facil- 
iate a structural transition from the pre-/3" to the (3" 
phase simply by providing the extra space needed for 
the (0,0,0) Mg atom to move 1/2 lattice vector in the 
z-direction. Including vacancies in the CE description is 
formally straight-forward, however, in practice it would 
increase the complexity in the CE and lead to other tech- 
nical complications. 
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VI. CONCLUSIONS 

We have investigated the formation of GP-zones in the 
Al-Mg-Si system using a CE of the total energy. DFT 
pseudo-potential calculations of 154 fcc-based pure, bi- 
nary and ternary structures were used to parameterize 
the energy by the SIM. Which cluster functions to in- 
clude was chosen according to the cross-validation crite- 
rion, with the intension to avoid over-fitting. The method 
successfully reproduces the total energy of 10 structures 
based on the suggested pre-/?" structure. Regarding the 
formation of GP-zones, the CE calculations shows that 
the MgiSii LIq structure suggested by Matsuda et al. 
is the most energetically stable phase. This is also con- 
firmed in direct DFT calculations. Our results are in 
fair agreement with the DFT calculations presented by 
Ravi and Wolverton,|3 where one should compare with 
their constrained GGA results for pre-/3" structures. Re- 



garding the formation of the pre-/?" and /?" structures, 
we outline two possible nucleation mechanisms. Either, 
the vibrational entropy {TS) term, or the presence of 
vacancies serve to stabilize the structures preceding the 
/?" phase, which then transforms to the meta-stable /?" 
phase. 
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Cluster, a Decorations, s 

units of (ao/2) 



Empty 














Point 
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i;2 






PI 


(000) 


(110) 




11; 12; 22 






P2 


(000) 


(200) 




11; 12; 22 






P7 


(000) 


(321) 




11; 12; 22 






Tl 


(000) 


(110) 


(101) 


111 


211 


222 






T2 


(000) 


(Oil) 


(Oil) 


HI 


211 


121 


221; 122; 


222 


T3 


(000) 


(110) 


(101) 


111 


211 


121 


221; 122; 


222 


T4 


(000) 


(121) 


(211) 


111 


211 


221 


222 




T5 


(000) 


(211) 


(112) 


111 


211 


121 


221; 122; 


222 


T6 


(000) 


(020) 


(112) 


111 


211 


221 


112; 212; 


222 


T7 


(000) 


(110) 


(112) 


HI 
112 


211 
212 


121 
122 


221; 
222 




Ql 


(000) 


(110) 


(101) (Oil) 


1111; 2111; 2 


211; 2221; 


2222 


Q2 


(000) 


(101) 


(Oil) (Oil) 


HI 


1; 21 


11; 2 


211; 1121; 


2121; 



2221; 1122; 2122; 2222 
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TABLE II: Interface energies between the Al phase and the 
Matsuda/pre-/3 phase. Energies in meV/A'^. We checked the 
negative value for the Al/pre-/3" (3 1 0) interface in a direct 
DFT calculation with result given in parenthesis. 

Al/pre-/3" 
(3 10) (2 3 0) (0 1) 



-3.57 (-4.31) 2.44 



6.44 



(10 0) 



Al/Matsuda 
(0 10) 



12.4 



5.9 



TABLE III: 



Structure 


pressure (meV/A^) 


Elastic energy 




(meV/A^) 




(mcV/B-phasc atom) 


4e (AlMgSia) 


-2.5777 




0.0488 


8cl (AlMg3Si4) 


8.9378 




0.5863 


Llo (MgSi) 


19.4610 




2.7798 


8cl (AlMg4Si3) 


28.1928 




5.8340 


4a (AlMgaSi) 


35.4267 




9.2119 


LI2 (AUMg) 


32.1312 




7.5778 



APPENDIX A: INPUT STRUCTURES 
APPENDIX B: PRE-/3 "-BASED STRUCTURES 



t-viV ^ 






^ / V / \ / 









FIG. 8: 32 atom structures included in the fitting procedure. 
Dashed rings represent Al atoms, and oversized/undersized 
rings represent Mg/Si atoms. New structures obtained by in- 
terchanging Mg and Si atoms were also included. Grey atoms 
are situated ao/2 above the white atoms. Half-circles indicate 
that two elements take up every second site in the vertical di- 
rection. 
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FIG. 9: 16 atom structures included in the fit. See further 
caption of Fig. |H| 
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FIG. 10: 12 atom structures included in the fit. See further 
caption of Fig. |H1 
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FIG. 11: 8 atom structures included in the fit. See further 
caption of Fig. |H1 
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FIG. 12: Some low-energy structures based on the pre-/3" 
structure. See further Fig. 2] 



